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Abstract 

We present calculations of heavy quark and quarkonium production at CERN LEP2 
in the /cT^-factorization QCD approach. Both direct and resolved photon contribution are 
taken into account. The conservative error analisys is performed. The unintegrated gluon 
distribution in the photon is taken from the full CCFM evolution equation. The tradi- 
tional color-singlet mechanism to describe non-perturbative transition of QQ-pair into a 
final quarkonium is used. Our analisys covers polarization properties of heavy quarkonia at 
moderate and large transverse momenta. We find that the total and differential open charm 
production cross sections are consistent with the recent experimental data taken by the L3, 
OPAL and ALEPH collaborations. At the same time the DELPHI data for the inclusive 
J/ ip production exceed our predictions but experimental uncertainties are too large to claim 
a significant inconsistency. The bottom production in photon-photon collisions at CERN 
LEP2 is hard to explain within the /cT-factorization formalism. 



1 Introduction 

Heavy quark and quarkonium production in photon-photon collisions at high energies is 
a subject of the intensive studies from both theoretical and experimental point [1-6]. From 
the theoretical point, heavy quark in 77 collisions can be produced via direct and resolved 
production mechanisms. In direct events, the two photons couple directly to a heavy quark 
pair. In resolved events, one photon ("single- resolved") or both photon ("double-resolved") 
fluctuate into a hadronic state and a gluon or a quark of this hadronic fluctuation takes 
part in the hard interaction. At LEP2 conditions the contribution of the double-resolved 
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events {gg QQ) is small [7], and charm and bottom quarks are produced mainly via 
direct (77 —>■ QQ) and single-resolved (751 — > QQ) processes. The direct contribution is 
not dependent on the quark and gluon content in the photon, whereas the single-resolved 
processes strongly depend on the photon's gluon density. Therefore detailed knowledge of 
gluon distributions in the photon is necessary for the theoretical description of such processes 
at modern (LEP2) and future (TESLA) colliders. 

Usually quark and gluon densities are described by the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parizi (DGLAP) evolution equation [8] where large logarithmic terms proportional 
to ln(//^) are taken into account. The cross sections can be rewritten in terms of process- 
dependent hard matrix elements convoluted with quark or gluon density functions. In this 
way the dominant contributions come from diagrams where parton emissions in initial state 
are strongly ordered in virtuality. This is called collinear factorization, as the strong ordering 
means that the virtuality of the parton entering the hard scattering matrix elements can be 
neglected compared to the large scale //^. However, at the high energies this hard scale is 
large compare to the Aqcd parameter but on the other hand fi^ is much less than the total 
energy ^/s (around 200 GeV for the LEP2 collider). Therefore in such case it was expected 
that the DGLAP evolution, which is only valid at large /x^, should break down. The situation 
is classified as "semihard" [9-12]. 

It is believed that at assymptotically large energies (or small x ~ A*^/s) the theoreti- 
cally correct description is given by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution 
equation [13] because here large terms proportional to ln(l/a;) are taken into account. Just 
as for DGLAP, in this way it is possible to factorize an observable into a convolution of 
process-dependent hard matrix elements with universal gluon distributions. But as the vir- 
tualities (and transverse momenta) of the propagating gluons are no longer ordered, the 
matrix elements have to be taken off-shell and the convolution made also over transverse 
momentum kj- with the unintegrated (/cj^-dependent) gluon distribution jF(a;, k|,). The un- 
intcgrated gluon distribution jF(x, k^) determines the probability to find a gluon carrying 
the longitudinal momentum fraction x and the transverse momentum ky. This generalized 
factorization is called /cr-factorization [10, 11]. 

The unintegrated gluon distribution is a subject of intensive studies [14]. Various ap- 
proaches to investigate this quantity have been proposed. One such approach, valid for both 
small and large x, have been developed by Ciafaloni, Catani, Fiorani and Marchesini, and 
is known as the CCFM model [15]. It introduces angular ordering of emissions to correctly 
treat gluon coherence effects. In the limit of asymptotic energies, it almost equivalent to 
BFKL [16-18], but also similar to the DGLAP evolution for large x and high /x^. The re- 
sulting unintegrated gluon distribution depends on two scales, the additional scale is a 
variable related to the maximum angle allowed in the emission and plays the role of the 
evolution scale /x^ in the collinear parton densities. We will use the following classification 
scheme [14]: ^(x, k|,) is used for pure BFKL- type unintegrated gluon distributions and 
A{x, k?j^, ii"^) stands for any other type having two scale involved. 

The CCFM evolution equation formulated for the proton has been solved numerically 
or analitically by different ways [19-21]. As it was shown [22-25], unintegrated gluon dis- 
tribution in the photon ^^(x, k|,, /x^) can be constructed by the similar method as in the 
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proton^. But situation is a bit different due to the pointlike component which reflects the 
sphtting of the photon into a quark-antiquark pair. Also in the photon there are no sum rules 
equivalent to those in the proton case that constrain the quark distributions. However, this 
difference is not significant because CCFM equation contains only gluon splitting g — > gg. 
For the first time the unintegrated gluon density in the photon was obtained [22] using a 
simplified solution of the CCFM equation in single loop approximation, when small-x effects 
can be neglected. It means that the CCFM evolution is reduced to the DGLAP one with 
the difference that the single loop evolution takes the gluon transverse momentum kx into 
account. Another simplified solution of the CCFM equation for a photon was proposed [23] 
using the Kimbcr-Martin-Ryskin (KMR) prescription [20]. In this way the /i^-dependence 
in the unintegrated gluon distribution enters only in last step of the evolution, and one-scale 
evolution equations can be used up to the last step. Both these methods give the similar 
results [23]. The phenomenological unintegrated gluon density, based on the Golec-Biernat 
and Wusthoff (GBW) saturation model [24] (extended to the large-x region), was proposed 
also [23]. Very recently the unintegrated gluon distribution in the photon A.y{x,h'^, n"^) was 
obtained [25] using the full CCFM evolution equation for the first time. It was shown that 
the full CCFM evolved effective (integrated over k|,) distribution is much higher than the 
usual DGLAP-based gluon density ai x < 10~^ region. 

In the previous studies [23, 25] the unintegrated gluon distributions in a photon (obtained 
from the single loop as well as full CCFM evolution equation) were applied to the calculation 
of the open charm and bottom production at LEP2. It was found that total cross section 
of the charm production is consistent with experimental data. In contrast, the theoretical 
predictions of the bottom total cross section underestimate data by factor 2 or 3. But we 
note that all these calculations reveal to the total cross sections of the open charm and 
bottom production only. In this paper we will study heavy flavored production at LEP2 
more detail using the full CCFM-evolved unintegrated gluon density [25]. We will investigate 
the total and differential heavy quark cross section (namely pseudo-rapidity and transverse 
momenta distributions of the D*-meson) and compare our theoretical results with the recent 
experimental data taken by the L3, OPAL and ALEPH collaborations at LEP2 [1-5]. 

Also we will study here the very intriguing problem connected with the inclusive J/ijj 
meson production at high energies. It is traces back to the early 1990s, when the CDF data 
on the J/ip and T hadroproduction cross section revealed a more than order of magnitude 
discrepancy with theoretical expectations. This fact has resulted to intensive theoretical 
investigations of such processes. In the so-called non-relativistic QCD (NRQCD) [28] there 
are additional (octet) transition mechanism from cc pair to the J/i/) mesons, where cc pair is 
produced in the color octet (CO) state and transforms into final color singlet state (CS) by 
help soft gluon radiation. The CO model describes well the heavy quarkonium production 
at Tevatron [29], although there are also some indications that it does not work well. For 
example, contributions from the octet mechanism to the J/ip photo- and leptoproduction 
processes at HERA are not well reproduce [30] experimental data. Also NRQCD is not 
predict J/ip polarization properties at HERA and Tevatron [30, 31]. At the same time usual 
CS model supplemented with fc^-factorization formalism gives fully correct description of 
the inclusive J/ip production at HERA [32] and Tevatron [33] including spin alignement of 

^See also [26], where we have used for unintegrated gluon density in a photon the prescription proposed 
by J. Bliimlein [27]. 
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the quarkonium. The CO contributions within the fc^-factorization approach also will not 
contradict experimental data if parameters of the non-perturbative matrix elements will be 
reduced [33, 34]. But this fact changes hierarchy of these matrix elements which are obtained 
within the NRQCD. 

Recently DELPHI collaboration has presented experimental data [35] on the inclusive 
J/ip production in 77 collisions at LEP2, which wait to be confronted with different theoret- 
ical predictions. The theoretical calculations [36] within the NRQCD formalism agree well 
with the DELPHI data. In contrast, the coUinear DGLAP-based leading order perturbative 
QCD calculations in the CS model significantly (by order of magnitude) underestimate [36] 
the data. The aim of this paper, in particular, is to investigate whether the inclusive J/ip 
production at LEP2 can be explained in the traditional CS model by using /cr-factorization 
and CCFM-based unintcgrated gluon density in a photon. 

The outline of this paper is following. In Section 2 we present the basic formahsm of the 
/cT-factorization approach with a brief review of calculation steps. In Section 3 we present 
the numerical results of our calculations. Finally, in Section 4, we give some conclusion. 
The compact analytic expressions for the off-shell matrix elements of the all subprocesses 
under consideration are given in Appendix. These formulas may be useful for the subsequent 
applications. 



2 Cross sections for heavy flavour production 

In 77 collisions heavy quark and quarkonium can be produced by one of the three mech- 
anisms: a direct production, a single-resolved and a double-resolved production processes. 
The direct contribution is governed by simple QED amplitudes which is independent on the 
gluon density in the photon. The double-resolved process gives a much smaller contribu- 
tion than the direct and single-resolved processes [7] and will not taken into account in this 
analysis. 

Let p^^^ and p^^^ be the momenta of the incoming photons and pi and p2 the momenta 
of the produced quarks. The single-resolved contribution to the 77 — >■ QQ process in the 
/cT-factorization approach has the following form: 

/dx — 
—A^{x, k^, ii'^)dk^—da{-fg* QQ), (1) 

where a{'yg* QQ) is the heavy quark production cross section via off-shell gluon having 
fraction x of a photon longitudinal momentum, non-zero transverse momentum kj^ (k|, = 
—k^^O) and azimuthal angle 0. The expression (1) can be easily rewritten as 

where \M.\'^{'yg* QQ) is the off-shell matrix element, s — {p^^ is the total cm. 

frame energy, 1/2 and 02 are the rapidity and azimuthal angle of the produced heavy quark 
having mass mq, 0.2 — m2T exp(|/2) / ^/s and m^x = '^q+pIt- To calculate the single-resolved 
contribution to the inclusive J/ip production the same formula (2) can be used where off- 
shell maxtrix element |A^P(75'* QQ) should be replaced by one which corresponds to the 
^yg* —>■ J/ip + g production process. 
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The available experimental data [1-5, 35] refer to heavy quark or quarkonium production 
in the e+e~ collisions also. In order to obtain these total cross sections, the 77 cross sections 
need to be weighted with the photon flux in the electron: 

d(T(e+e~ e^e~QQ + X) ^ j f^/e{xi)dxi j f^/e{x2)dx2 (^(7(77 QQ), (3) 

where we use the Weizacker- Williams approximation for the bremsstrahlung photon distri- 
bution from an electron: 

with Ql,^ = mlxy{l - xf and Ql,^ = Ql^ + {Eef{l - x). Here x = E^/E„ E = E, ^ 
\/s/2, s is the total energy in the e+e" collision and ^ ~ 30 mrad is the angular cut that 
ensures the photon is quasi-real. 

The multidimensional integration in (2) and (3) has been performed by means of the 
Monte Carlo technique, using the routine VEGAS [37]. The full C-|— I- code is available from 
the authors on request^. For reader's convenience, we collect the analytic expressions for the 
off-shell matrix elements \A4\'^{'~fg* QQ) and \A4\'^{jg* J/ipg) in Appendix, including, 
in particular, the relevant formulas for helicity zero J/ip production state. 



3 Numerical results 

First of all, we can perform integration of the unintegrated gluon distribution 
Ay(x, k^, fji^) [25] over gluon transverse momenta to obtain the effective gluon density in 
the photon: 

xg^{x,ix^) ~ j ^^(x, k^, //.^)(ik^. (5) 


The effective density xg^{x,n'^) can be compared with the experimental data [38, 39] taken 
by HI collaboration at HERA. As seen in Figures 1 and 2, this gluon distribution agrees 
well to the existing data extracted from the hard dijet (mean p'^ = 38 GeV^ [38] and p'^ — 
74GeV^ [39]). In contrast, KMR construction of unintegrated gluon density in the photon 
tends to underestimate the HERA data at a; < 10^^ [23]. 

Being sure that the full CCFM-evolved unintegrated gluon density Ay{x,k'^, n'^) repro- 
duces well the experimental data for the xg-y{x, /i^), we now are in a position to present our 
numerical results. We describe first our theoretical input and the kinematical conditions. 
The cross sections for heavy quark and quarkonium production depend on the heavy quark 
mass and the energy scale /x^. Since there are no free quarks due to confinement effect, 
their masses cannot be measured directly and should be defined from hadron properties. In 
our analysis we have examined the following choice: rric — 1.4 ± 0.1 GeV for charm and 
m^j = 1.75 ± 0.25 GeV for l)ottom quark masses. Such variation of the quark masses gives 

^ lipatov@t heory. sinp . msu 
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the largest uncertainties in comparison with scale variation'^ and therefore can be used as 
an estimate of the total theoretical uncertainties. Then, we will apply standard expression 
A*^ — + Pt for both renormalization and factorization scales. Here Pt is the transverse 
momentum of the heavy quark in the center-of-mass frame. We use LO formula for the 
strong coupling constant with Uf = 4 active quark flavours and Aqcd — 200 MeV, 

such that as(M|) = 0.1232. But Aqcd = 340 MeV [7] was tested also. 

Figure 3 confronts the total cross section 17(77 ~^ cc + X) calculated as a function of the 
photon-photon total energy W^^ with experimental data [2] taken by the L3 collaboration 
in the interval 5 < Wjj < 70 GeV. Sohd line represent the calculations with the charm 
mass rric = 1.4 GeV, whereas upper and lower dashed lines correspond to the rric = 1.3 
GeV and rric = 1.5 GcV respectively. It is clear that /cx-factorization reproduces well both 
the energy dependence and the normalization. One can see that sensitivity of the results to 
the variations of the charm mass is rather large: shifting the mass down to mc = 1.3 GeV 
changes the estimated cross section by 15 — 20% at W^j ~ 60 GeV. But in general all three 
curves lie within the experimental uncertainties. 

Experimental data for the total cross section a{e^e~ e^e~cc + X) come from the 
three LEP collaborations L3 [1, 3], OPAL [4] and ALEPH [5]. In Figure 4 we show our 
predictions in comparison with data. All curves here are the same as in Figure 3. One can 
see that our calculations describe the experimental data well again. The variation of the 
quark masses 1.3 < rric < 1-5 GeV gives the theoretical uncertainties approximately 15% in 
absolute normalization. 

The available experimental data were obtained for the D* meson production also. Two 
differential cross section are determined: the first one as a function of the transverse D* 
momentum pt, and the second as a function of pseudo-rapidity \r]\. In our calculation we 
convert charmed quark into D* meson using the Peterson fragmentation function [41]. The 
default set for the fragmentation parameter Cc and the fraction /(c ^ D*) is Cc = 0.06 
and /(c D*) = 0.26. Other values for the parameter ec in the NLO perturbative QCD 
calculations are often used also, namely ec = 0.116 [42] in the massless scheme and ec = 
0.031 [43] in the massive one. In the case of the massless calculation, this parameter was 
determined via a NLO fit to LEPl data on D* production in e~^e~ annihilation measured 
by OPAL collaboration [44]. To investigate the sensitivity of the our numerical results to Cc 
parameter wc have repeat our calculations using ec = 0.031. 

The recent L3 data [3] refer to the kinematic region defined by 1 < < 12 GeV and 
\r)\ < 1.4 with averaged total e+e" energy ^/s = 193 GeV (183 < i/i < 209 GeV). The 
OPAL data [4] were obtained in the region 2 < < 12 GeV, \n\ < 1.5 and averaged 
over 183 < y/s < 189 GeV. The more recent ALEPH data [5] refer to the same kinematic 
region but averaged over 183 < ^/s < 209 GcV. We compare these three data sets with our 
calculation at y/s = 193 GeV. The different values of y/s are not expected to change the 
cross section more than the corresponding experimental errors. We have checked directly 
that shifting from 189 to 193 GeV increase the calculated cross sections by about one 
percent only. 

The transverse momenta distributions of the D* meson for different pseudo-rapidity re- 
gion in comparison to experimental data shown in Figure 5 and 6. The solid and both 

"^It was shown [25] that Monte Carlo generator CASCADE [19, 40] predicts the very similar results for 
charm total cross section with both scales /U^ = rrig -|- p|. and /i^ = ^tuq. 



6 



dashed curves here are the same as in Figure 3 (calculated with the default value ec = 0.06), 
dash-dotted curves represent results obtained using eg = 0.031 and rric = 1.4 GeV. The 
overall agreement between the our predictions and the data is good although the ALEPH 
data points in medium and large pt range lie slightly above the theoretical curves. Shifting 
the Sc default value down to = 0.031 results to a bit broadening of the px spectra, which is 
insufficient to describe the data. The effects come from changing of the charm mass present 
only at low pt- the predicted cross section da/dpT with rric = 1.3 GeV is 10 — 15% above the 
one calculated with rric — 1 A GeV at pt ~ 1 GeV, whereas sohd and both dashed theoreti- 
cal curves practically coincide at medium and large pt- The similar effect was found in the 
NLO perturbative calculation [45], where the difference between the massive and massless 
approach arises at low pt only. 

The D* pseudo-rapidity distributions compared with the experimental data in different 
Pt range are shown in Figure 7 and 8. All curves here are the same as in Figure 5. Our 
calculations agree well with measured differential cross sections but slightly underestimate 
the OPAL data. However, setting = 0.031 increases the absolute normalization of the 
pseudo-rapidity distribution by approximately 10%, and agreement with OPAL data becomes 
better. Again, one can see that the significant mass dependence has place at low transverse 
momenta only: the difference between the theoretical curves calculated at 2 < < 12 GeV 
and plotted in Figure 8 is much smaller than difference between the results presented in 
Figure 7 obtained at 1 < p^ < 12 GeV. 

We conclude from Figs. 3 — 8, that our calculations agree well with charm data at 
LEP2. In contrast with charm case, the open bottom production in 77 collisions is clear 
underestimated by /cx-factorization. Figure 9 shows our prediction for the open bottom 
cross section compared to the L3 [1] and OPAL [4] experimental data. Using low but still 
reasonable 6-quark mass = 4.5 GeV, we obtain a{e^e^ e^e^bb + X) = 2.94 pb at 
^/s = 200 GeV. The very similar value a = 2.7 pb was obtained [23] within the GBW 
saturation model adopted for the photon. The Monte Carlo generator CASCADE predicts 
a — 4.9 pb [25] where normalization factor n = 1.7 has been applied for the resolved 
contributions. The calculation [23] based on the KMR prescription for unintegrated gluon 
density gives a lower cross section a = 1.9 pb. At the same time the prediction of the 
massive NLO perturbative QCD calculation [43] is 3.88 pb for = 4.5 GeV and 2.34 pb 
for nib = 5.2 GeV respectively. All these results are significantly (about factor 2 or 3) lower 
than experimental data. 

Such disagreement between theory and data for bottom production at LEP2 is surprising 
and needs an explanation. It is known that the similar difference between theory and data 
was claimed for inclusive bottom hadroproduction at Tevatron. Recent analysis indicates 
that the overall description of the these data can be improved [46] by adopting the non- 
perturbative fragmentation function of the 6-quark into B meson: an appropriate treatment 
of the 6-quark fragmentation properties considerably reduces the disagreement between mea- 
sured bottom cross section and the corresponding NLO calculations. It would be interesting 
to find out whether the similar explanation is also true for the L3 and OPAL experimental 
data. 

After we have studied open heavy quark production at LEP2, we will investigate pro- 
duction of the heavy quarkonum in 77 interactions. As it was already mentioned above, 
non-relativistic QCD gives a good description [36] of the recent DELPHI data [35] on inclu- 
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sive J/ip production at LEP2. We will examine whether the DELPHI data can be expained 
within the CS model using /c^-factorization approach and CCFM-based unintegrated gluon 
density in a photon. Again, only direct and single-resolved contributions are taken into 
account. 

Now we change the default set of parameters which were used in the case of open charm 
calculations. Since we neglect the relative momentum of the c-quarks (which form a, J/ip 
meson) the charm mass should be taken rric = Therefore as default choice in following 

we will use rric — 1.55 GeV. On the other hand there are many examples when smaller value 
rric — 1.4 GeV in the calculation of J/t/j production is used [30, 32, 36]. In our analysis we 
will apply this value as extremal choice to investigate the theoretical uncertainties of the 
calculations. 

The DELPHI data [35] refer to the kinematic region defined by — 2 < < 2 with total 
e+e~ energy y/s — 197 GeV, where is the J/ip rapidity. However these data were obtained 
starting from very low transverse momenta p^j^ ~ 0.2 GeV^. We note that /cr-factorization 
as well as usual coUinear factorization theorem does not work well for such p^T values, and 
our calculations should be compared with experimental data at approximately p^T > 1 GeV 
only. In Figure 10 we confront our theoretical predictions with the measured differential 
cross section da/dp^j,. Solid line corresponds the default set of parameters and lies below 
the data by a factor about 2 or 3. This dicrepancy is not catastrofic, because some reasonable 
variations in rric and Aqcd, namely 200 < Aqcd < 340 MeV and 1.4 < rric < 1-55 GeV, 
change the estimated cross section by a factor of 3 (dashed line in Figure 9). Thus the 
visible disagreement is eliminated. However, we do not interpret this as a strong indication 
of consistency between data and theory, but rather as a concequence of a wide uncertainty 
band. Better future experimental studies are crucial to determine whether the results of our 
calculations contradict DELPHI data points. 

The main difference between /c^-factorization and other approaches connects with po- 
larization properties of the final particles because the initial off-shell gluons do promptly 
manifest in the J/ip spin alignement [32-34]. Only a very small fraction of J/ip mesons can 
be produced in the helicity zero state (longitudinal polarization) by massless bosons. This 
property is totally determined by the subprocess matrix clement structure. The degree of 
spin alignement can be measured experimentally since the different polarization states of 
result in significantly different angular distributions of the ^ 1^^ decay leptons: 

"^W-f^')~l + acos-9. (6) 
acosb* 

where a = 1 for transverse and a = —1 for longitudinal polarizations, respectively. Here 9 
is the angle between the lepton and J/V^ directions, measured in J/'^ meson rest frame. We 
calculate the p|.-dependence of the spin alignement parameter as 



1 - 3e(p^T) 
1 + ^(p^t) 



with C(p^r) being the fraction of longitudinally polarized J/ if) mesons. The results of our 
calculations are shown in Figure 11. Solid line represents the A;T-factorization predictions 
and dashed line corresponds to the coUinear leading-order pQCD ones with GRV (LO) 
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gluon density [47] in a photon. One can see that fraction of longitudinally polarized mesons 
increases with within /cx-factorization approach. This fact is in clear contrast with usual 
coUinear parton model result. The A;T-factorization calculations made for the inelastic J/'^ 
production at HERA [32, 34] and Tevatron [33] posses the same behavior of the a parameter, 
whereas collinear parton model and NRQCD predict the strong transverse polarization at 
moderate and large range. We point out that our predictions for the J j^il) polarization are 
stable with respect to variation of the model parameters, such as charmed quark mass and 
factorization scale. In fact, there is no dependence on the strong coupling constant which is 
cancels out. At the same time the DELPHI fit [35] gives a — 0.7±1.3 for p^-^ > 1 GeV^, which 
has, however, huge experimental uncertainties. Since account of the octet contributions does 
not change predictions of the fcx-factorization approach for the spin alignement parameter 
a [33], the future extensive experimental study of such processes will be direct probe of the 
gluon virtuality. 

4 Conclusions 

We have investigated a heavy flavour production in photon-photon interactions at high 
energies within the framework of /c^^-factorization, using unintegrated gluon distribution 
obtained from the full CCFM evolution equation for a photon. We calculate total and 
differential cross sections of the open charm and bottom production including D* meson 
transverse momenta and pseudo-rapidity distributions. Also we have studied inclusive Jjil) 
production at LEP2 using color-singlet model supplemented with /sT-factorization. 

We take into account both the direct and single-resolved contribution and investigate the 
sensitivity of the our results to the different parameters, such as heavy quark mass, charm 
fragmentation and Aqcd parameter. There are, of course, also some uncertainties due to 
the renormalization and factorization scales. However, these effects would not to be large 
enough to change the conclusions presented here, and were not taken into account in our 
analisys. 

The results of calculations with default parameter set agree well with open charm pro- 
duction data taken by the L3, OPAL and ALEPH collaborations at LEP2. In contrast, 
bottom production cross section is clearly underestimated by a factor about 3. A potential 
explanation of this fact may be, perhaps, connected with the more accurate treatment of the 
6-quark fragmentation function. Our prediction for the inclusive J/V^ production slightly un- 
derestimate the DELPHI data. However, a strong inconsistensy cannot be claimed, because 
of large experimental errors and theoretical uncertainties. Therefore more precise future ex- 
periments, espesially polarized quarkonium production, are necessary to know whether our 
predictions contradict the experimental data. 

In conclusion, we point out that at CERN LEP2 collider (as well as at HERA and 
Tevatron) the difference between predictions of the collinear and /c^-factorization approaches 
is clearly visible in polarized heavy quarkonium production. It comes directly from initial 
gluon off-shellness. The experimental study of such processes should be additional test of 
non-coUinear parton evolution. 
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6 Appendix 

Here we present compact analytic expressions for the off-shell matrix elements which 
appear in (2). In the following, s, t, u are usual Mandelstam variables for 2 — > 2 process and 
cq is the fractional electric charge of heavy quark Q. 

We start from photon-gluon fusion 'y{q)g*{k) — > Q{pi)Q{p2) subprocess, where initial off- 
shell gluon has non-zero virtuality k^. The corresponding squared matrix element summed 
over final polarization states and averaged over initial ones read 

\Mn^9* - QQ) - _ FQ^(ki), (Ai) 

{t — m)[u — m) 

where m is the heavy quark mass, and 

FQQ(k|) = 6m^ - (2k^ + 2{i + + 3P + Sm^ + Uiu)n/+ 

{2{i + u)]4- + Sivkl + + u^ + Ttu^ + 7Pu)m'^- 

i^2]4' + 2{i + u)\4 + P + u^)- {A.2) 

We note that matrix element of the direct contribution 77 — > QQ may be easily obtained 
from (A.l) and (A. 2) in the limit k|, — > 0, if we replace normalization factor {Anyeqaas 
by the {A7T)'^a'^eQNc (where N^. is the number of colors) and average (A. 2) over transverse 
momentum vector k^. 

Now we are in a position to present our formulas for the 'j{q)g*{k) — > J/ip{p^)g{pg) 
subprocess. In the color-singlet model the production of J/ip meson is considered as pro- 
duction of a quark-antiquark system in the color-singlet state with orbital momentum L = 1 
and spin momentum S = 1. The squared off-mass shell matrix element summed over final 
polarization states and averaged over initial ones can be written as 

, , . 64e^(47r)W|^(0)P , , , 

3m^(i - m^y{u - - k^)^(i + u + k^yk^ 

where |^(0)p = 0.0876 GeV^ is the J/ip wave function at the origin, is the J/i/j meson 
mass, k^ is the virtuality of the initial gluon, and function F^(ky) is given by 

F^(k|) = k|(k^(m5, - i){ml - i - u) - ml{P + iu + u^- 

ml{i + u)f + k^(3m5, + iu{i + - 3mJ(2£ + -u) + mj(3? + 2iu - iy))+ 
kl{2ml + mti{i - u) - P{i + uf - 2m^(2t + u) + mli{2p + 5iu + 5u^))) + 
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2k^(k^(mJ - t)(mj -i-u)+ l4{3ml + i{i + 11)^- mj(5i + 3u)+ 

+ 2iu + 3?i^)))(-|p^T||kT| cos^a) + 2ml{ml{-ml + i + uf+ 

k|(mj + P + 2£n - - 2mJ(t + M)))p2^^k| cos^ 02. (A.4) 

Here p^y is the J /ip transverse momentum, 02 is the azimutal angle of the incoming virtual 
gluon having virtuality k^. To study polarized J/ if) production we introduce the four- 
vector of the longitudinal polarization e^^. In the frame where the z axis is oriented along 
the quarkonium momentum vector, = (£"^,0,0, |pv>|), this polarization vector is e^^ = 
(Ip^I, 0, 0, E^)/m^. The squared off-shell matrix element read 

where function F^^L(k|,) is defined as 

F^A^^) = -4A2A3k^mJ + 2A^k^mJ - 4A2A3k| 

8A2 Agk^mjt - 4A^k^mJt + 8A2 Agklmjt + 4A^k^mJt- 
4A2A3k^P + 2A^k^P - 4A2A3k^m5,? - SA^k^mJ? + 2A^mJP+ 
4A2k2 p _ 4A2m2 p + 2A2£4 + ^^^\m\u- 
4A2A3k|,m5,it 2A^k|,m5,rt 4A^k|mJit 2A^k|,mJ,it- 
4A^k^£M 4A2A3k^t'u - 2A^k^tM - AAlklmltu - 4A2A3k|mJtM- 
4A^k|mJ£M + 2A^mJ,£M + 4A2A3k^pM + 2A^k^pM - 6A^mJ?M+ 
4A3F'a k^M^ - AAlklmlu^ + 4A2 A3kymJ'a^ - 2Al\4'mlu^ + k^mJuV 

4A2A3k|iV - 2A3ky£^i^ - 2k^iV - 2Almliu'^ - 2k^mliu^+ 
2Alt^u^ + 2\4Pu^ - 2k^mJ# + 2kliu^ + k^u"^ + SAjk^{-ml + i + uf+ 
4Aik|(A2(k^(2mJ - 2i - u) + m\{m\ - 2i - u)u + k^(2mj + -u^- 
2ml{i + u))) + Asi-ml + k^{-ml + t) + 2i{i + uY+ 
ml{3i + u) - mli{4i+3u) + k|(-(mj(2£ + + t{2t + 3u))))+ 
4(-(A2A3(k2 + ml){ml - £)(mj - t - u)) + k|(k^ - i>V 
4A2k|(-mJ + i + uf + Ai(2A2k|(2mJ + k^(2mj - 2t - u)- 
2ml{i + u)+ u{i + u)) - A3(m^ + k^(mj - t + «)- 
2i(£ + - mj,(3i + 2?i) + mj(4p + 5i?i + u^)+ 
k|(-2P - tii + ii^ + m|(2£ + 'ji)))))(-|p^T||kT| cos02)+ 
4(2AiA2(k| + mj)(mj - t - u) + 
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+ 2Al{-ml + i + uf)plT'^lcos^ 02, (^.6) 
and the following notation has been used: 



4Jp2^ + s(ai-A)V4 



A3 = =|p^T||Kr| COS02, 



"1 = \/ — exp(7/^), /3i = \ ^ exp(-7/^), 
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Figure 1: The effective gluon distribution xg^{x,n'^) in the photon as a function of log^g^ 
at fji^ — 38GeV^. The experimental data [38] taken by HI collaboration from hard dijets 
production at HERA. 
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Figure 2: The effective gluon distribution xg^{x,fj,'^) in the photon as a function of x at 
//^ = 74GeV^. The experimental data [39] taken by HI collaboration from hard dijets 
production at HERA. 
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Figure 3: The open charm total cross section o-(77 cc + X) as a function of W^j at 
^/s = 189 — 202 GeV. The sohd hne corresponds to the charm mass rric = 1.4 GeV, upper 
and lower dashed lines correspond to the rric = 1.3 GeV and rric = 1.5 GeV respectively. 
The experimental data [2] are taken by the L3 collaboration. 
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Figure 4: The open charm total cross section ^(e+e — > e+e cc + X) as a function of the 
e'^e~ centre-of-mass energy y/s. All curves are the same as Figure 3. The experimental data 
are from L3 [1, 3], OPAL [4] and ALEPH [5] collaborations. 
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Figure 5: The differential cross section da/dpT for inclusive D* production at \ri\ < 1.4. 
The solid and both dashed curves here are the same as in Figure 3 (calculated with the 
default value eg = 0.06), dash-dotted curve represents results obtained using Cc = 0.031 and 
mc = 1.4 GeV. The experimental data are from L3 [3]. 
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Figure 6: The differential cross section da/dpr for inclusive D* production at I77I < 1.5. All 
curves here are the same as in Figure 5. The experimental data are from OPAL [4] and 
ALEPH [5]. 
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Figure 7: The differential D* cross section da/d\ri\ for the process e'^e~ — > e^e~D* -\- X va. 
the 1 < pt < 12 GeV range. All curves here are the same as in Figure 5. The experimental 
data are from L3 [3] . 
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Figure 8: The differential D* cross section da/d\r]\ for the process e'^e" — > e+e~D* + X in 
the 2 < pt < 12 GeV range. All curves here are the same as in Figure 5. The experimental 
data are from OPAL [4] and ALEPH [5] collaborations. 



22 




140 160 180 200 220 

Vs (GeV) 

Figure 9: The bottom total cross section a{e'^e~ e^e~hh + X) as a function of the 
e^e~ centre-of-mass energy ^/s. The sohd hne corresponds to the bottom mass mf, = 4.75 
GeV, upper and lower dashed lines correspond to the = 4.5 GeV and rub = 5.0 GeV 
respectively. The experimental data are from L3 [1] and OPAL [4]. 
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Figure 10: The cross section da/dp^j, of e+e~ e'^e~J/ijj + X measured by DELPHI 
collaboration [35] at 2 < < 2 as a function of p^prp is compared with the /cr-factorization 
calculations in the CS model. The solid line corresponds to the rric = 1.55 GeV and Aqcd = 
200 MeV, dashed hne corresponds to the rric^lA GeV and Aqcd = 340 MeV. 
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Figure 11: The p^rp dependence of the spin ahgnement parameter a for the inclusive J/'^ 
production. The sohd hne corresponds to the /c^-factorization predictions, dashed hne cor- 
responds to the colhnear leading order pQCD calculations with the GRV (LO) gluon density 
in a photon. 
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